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ABSTRACT 

We present Gemini near-infrared spectroscopic observations of six luminous quasars 
at z = 5.8 ~ 6.3. Five of them were observed using Gemini-South/GNIRS, which 
provides a simultaneous wavelength coverage of 0.9-2.5 /im in cross dispersion mode. 
The other source was observed in K band with Gemini-North/NIRI. We calculate line 
strengths for all detected emission lines and use their ratios to estimate gas metallicity 
in the broad- line regions of the quasars. The metallicity is found to be supersolar with 
a typical value of ~ 4 Z©, and a comparison with low-redshift observations shows no 
strong evolution in metallicity up to z ~ 6. The Feii/Mgll ratio of the quasars is 
4.9 ± 1.4, consistent with low-redshift measurements. We estimate central BH masses 
of 10 9 to 10 10 Mq and Eddington luminosity ratios of order unity. We identify two 
Mgn AA2796,2803 absorbers with rest equivalent width W A2796 > 1 A at 2.2 < z < 3 
and three Mgn absorbers with VF^ 2796 > 1.5 A at z > 3 in the spectra, with the two 
most distant absorbers at z = 4.8668 and 4.8823, respectively. The redshift number 
densities (dN/dz) of Mgn absorbers with W^ 2796 > 1.5 A are consistent with no cosmic 
evolution up to z > 4. 
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INTRODUCTION 



In recent years we have wi tnessed the discoveries of quasars at z ~ 6; at this epoch the universe 



was less than one Gyr old (e.g. iFan et al.112003 . 120061 ; ICool et al. 



2006 



McGreer et alJl200d ). These 



high-redshift quasars are among the most luminous objects known and provide direct probes of 
the early universe when the first generations of galaxies and quasars formed. Quasar activity 
and the formation processes of galaxies and supermassiye black holes (BHs) are closely correlate d 



(e.g. iKauffmann &i Haehnelt 



2000 



Wvithe & Loeb 



2003 



Hopkins et al J 120051 : ICroton et al J 12006 ) 



so these z ~ 6 objects are essential for studying the accretion history of BHs, galaxy formation, 
and chemical evolution at very early epochs. High-redshift luminous quasars harbor BHs with 
masses of 1 9 -10 1 M M and em it near the Eddington limit (e.g. iBarth et al.ll2003l ; IVestergaardll2004l ; 
Jiang et al.ll2006l ; iKurk et al.ll2007l ). revealing the rapid growth of central BHs. Their emission line 



ratio s show solar or super s olar metallicity, as is found in low-redshif t quasars (e.g. IBarth et al 



20031 ; iDietrich et al.1 l2003al ; iFreudling et al.l 120031 ; iMaiolino et al.l 120031 ) , indicating that there was 
vigorous star formation and element enrichment in the first Gigayear. 

Understanding quasars requires observations from X-ray to radio wavelengths; each spectral 
region is due to different physical mechanisms and probes different regions of the active nucleus. In 
particular, the rest-frame UV spectrum, from the Lya emission line to the Fe n bump at 2000-3000 
A, contains strong diagnostic emission lines and provides key information on the physical conditions 
and emission mechanisms of the broad-line region (BLR). The bulk motions of the BLR can be used 
to determine the mass of the central BH, while chemical abundances in the BLR are important in 
understanding the history of star formation in the host galaxy. Photoionization models show that 
various emission-line ratios provide re liable estimates of metallicity in the BLRs of quasars (e.g. 
Hamann et al.ll2002l : lNagao et alJbOQfl j. Strong permitted line ratios, such as N v A1240/C iv A1549 
(hereafter Nv/Civ), Nv A1240/Heii A1640 (hereafter Nv/Heii), and F ell/Mgil A2800 (hereafter 
Fell/Mg n), are relatively easy to measure even in z ~ 6 quasars (e.g. 



199 



Hamann fc Ferland 
19931. 



1993 



most of the 



19991 ) . Fe abundance is of particular interest (e.g. lHamann Sz Ferlandl 
Fe in the solar neighborhood is generated by Type la supernovae (SNe la), whose precursors are 
believed to be long-lived intermediate-mass stars in close binaries, so appreciable Fe enrichmen t 
can only happen on a timescale of one Gyr after the initial starburst (e.g. iGreggio fc Renzinil ll983). 
On the contrary, the production of a elements such as Mg and O is dominated by SNe of types II, 
lb and Ic, which explode very soon after the initial burst. Therefore, the Fe/a ratio is expected to 
be a strong function of age in young systems, and thus it puts a useful constraint on the age and 
chemical enrichment of the gas in the quasar environment. 

At low redshift, elemental abundances measured from strong emission and intrinsic absorption 
lines have revealed that quasar environments have solar or supersolar metallicity. At z ~ 6, 



NIR range. NIR spectroscopy has been o b tained for a few z 



Pentericci et al.1 hood: barth et al.l I2OO3I : IFreudling et al.l I2OO3I : IMaiolino etal 



2003 



Iwamuro et al 



quasars (e.g. 


Goodrich et al. 


2001; 


Maiolino et al. 


2003 


> 


Stern et al. 



2004J). These observations only provide NIR spectra with low signal-to-noise 
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ratios (SNRs) or partial wavelength cov erage. Even basic p roperties such as the average continuum 
slope are not well measured. Recently, iKurk et al.l (|2007l ) obtained NIR spectra with high SNRs 
for a sample of five z ~ 6 quasars with ISAAC on VLT. They used the ISAAC data to study the 
Feli/Mgil ratios, BH masses, and the Stromgren spheres around quasars. In this paper we present 
Gemini NIR spectroscopy of six z ~ 6 quas ars. All these quasars were discovered by the Sloan 



Digital Sky Survey (SPSS; lYork et alJl200(J ) and were selected from Fan et al. 2000, 2001, and 
2004 (hereafter quasa r discovery papers). T hey have extensiv e follow-up observations in X-ray (e.g. 
Shemmer et al.l 120061'). Svitzer ban ds (e.g. I Jiang et al.l 120061 ) . millimeter /submillimeter and radio 



wavelengths (e.g. IWang et al.ll2007l ). providing a fundamental dataset to study luminous quasars at 



z ~ 6 in the context of the growth of early BHs and their relation to galaxy formation. The Gemini 
data have high SNRs and excellent NIR wavelength coverage. They are used to measure the metal 
abundances in high-r edshift environments and estimate central BH masses using empirical mass 
scaling relations (e.g. iMcLure fc Dunlopl 12004 ; IVestergaard fe Peterson! 120061 ) . 



We also use these high-redshift spectra to study the evolution of strong intergalactic Mg II ab- 
sorption at z < 6. Quasar absorption lines are an unbiased tracer of intervening gaseous material 
along the lines of sight to quasars. Strong Mgll absorption systems are usually d i rectly associ 



1998 



Churchill et al. 



2000 



Steidel et al 



ated with galaxy halos or disks (e.g. ICharlton Churchill 
2002). The statistical properties of the Mgn AA2796,2803 doublet absorbers at low redshift 



(z < 2.2) have been well studied (e.g. Steidel & Sargent 1992; Nestor et al. 2005; Prochter et 
al. 2006; see a summary of Mgn absorption surveys by Churchill et al. 2005). At higher redshift, 
the Mgll doublet is shifted to the NIR wavelength range, where detectors are much less efficient 
than are detectors at the optical region. The study of high-redshift Mgn absorbers is therefore 
limited by the small number of quasars with good NIR spectra available. Our Gemini/GNIRS 
spectra with high SNRs have excellent NIR wavelength coverage, so they are efficient to probe the 
Mgn absorption systems at z < 6. 

The structure of the paper is as follows. In § 2 we describe our Gemini observations and 
data reduction of the six quasars. In § 3 we conduct a detailed spectral analysis, including the 
measurement of emission- line strengths. Chemical abundances are calculated and central BH masses 
are estimated in § 4. Mgn absorption systems in the quasar spectra are identified and analyzed 
in § 5. We give a brief summary in § 6. In this paper we use A (Ao) to denote observed-frame 
(rest-frame) wavelen gth. We use a A-do minated flat cosmology with Ho = 70 km s -1 Mpc -1 , fl m 
= 0.3 and tt A = 0.7 (jSpereel et aHl2007h . 



2. OBSERVATIONS AND DATA REDUCTION 

We obtained NIR spectra for the six luminous z ~ 6 quasars using Gemini-South/GNIRS and 
Gemini-North/NIRI. Five of them were chosen to be easily reached by Gemini-South, and were 
observed using GNIRS in February and March 2006. The observations on GNIRS were carried out 
in cross-dispersion mode, giving a simultaneous wavelength coverage of 0.9-2.5 jivci with excellent 
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transmission. At z ~ 6, strong diagnostic emission lines such as Civ and Mgn are redshifted 
to this wavelength range, so GNIRS provides us an efficient way to collect NIR information for 
high-redshift quasars. We used the short camera on GNIRS with a pixel scale of 0.15"/pixel. The 
slit length in cross-dispersion mode is 6.1" and we used a slit width of 0.675" to optimize the 
combination of resolution and light gain. The resolving power is about 800 and the dispersion 
varies from ~3 A per pixel at 0.9 /im to ~6 A per pixel at 2.5 /im. This is sufficient to separate 
the Mgn AA2796,2803 doublet absorption lines. We used the standard ABBA nodding sequence 
between exposures. The exposure time at each nod position was five minutes and the distance 
between the two positions was 3". Before or after the exposure of each quasar, a nearby A or F 
spectroscopic standard star was observed for flux calibration and to remove telluric atmosphere 
absorption. The log of observations is given in Table 1, where redshifts, zab (AB magnitudes) and 
J (Vega-based magnitudes) are taken from the quasar discovery papers. 

The GNIRS data were reduced using the Gemini package within IRAE^l Briefly, for each 
quasar the NIR spectroscopic data were sky-subtracted for each pair of images taken at the two 
nod positions. Then the images were flat-fielded. In cross-dispersion mode, GNIRS uses orders 3-8 
to cover 0.9-2.5 /im. A single flat has non-uniform illumination between orders, so three sets of 
flats with different exposure times were taken each night. The final flat is extracted from the three 
flats so that each order has good illumination but without saturation. After flat-fielding, distortion 
correction and wavelength calibration were applied, the frames were combined, and one-dimensional 
spectra were extracted. Then the spectra were flux calibrated and corrected for telluric atmosphere 
absorption using the spectra of the spectroscopic standard stars. Finally the spectra in different 
orders were scaled and connected to a single spectrum. Adjacent orders cover a short common 
wavelength range, which enables the individual orders to be scaled to form a continuous, smoothly 
connected spectrum. This arbitrarily scaled spectrum was then calibrated to the J broad-band 
magnitude listed in Table 1 and thereby placed on an absolute flux scale with an uncertainty of 
5~10%. The typical SNR at 1.2 /im reaches ~ 8 per pixel. 

Figure 1 shows the spectra of the five quasars in the observed frame. The spectra at A < 1.0 
jj,m are taken from the quasar discovery papers. The spectra have been smoothed by three pixels. 
The bottom panel shows the sky transparency at an airmass of 1.0 and water vapor of 1.6 mm. 
There are strong telluric absorption bands around 1.4 and 1.9 /mi, so we do not show the spectra 
in these ranges. In Figure 1 we also show a spectrum of a standard star HIP 54027 calibrated using 
the same procedure as we do for our science objects. For comparison, a spectrum (in red) of a black 
body with the effective temperature of HIP 54027 is given on top of the spectrum of HIP 54027. 
The good agreement between the two spectra indicate that spectra in different dispersion orders 
are properly scaled and connected. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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We obtained a K-band spectrum for SDSS J1623+31lfl using Gemini-North/NIRI in August 
and September 2004 and August 2005. The observations were made with a f/6 camera in spectro- 
scopic mode. The long-slit has a length of 110" and we chose to use a slit width of 0.75", which 
provides a resolving power of ~ 500. The observing strategy and data reduction are similar to those 
used for the GNIRS data. The standard ABBA mode was used with a nodding distance of 20" and 
an exposure time of five minutes at each nod position. A nearby A or F spectroscopic standard star 
was observed before or after the target observation. The NIRI data were also reduced using the 
Gemini package within IRAF. After sky-subtraction, flat-fielding and wavelength calibration, the 
frames were combined. Then one-dimensional spectra were extracted from the combined images, 
and were flux calibrated using the spectra of standard stars. The observing log is given in Table 1 
and the spectrum is shown in Figure 2. Since we do not have K-b&nd photometry for this object, 
we cannot apply an absolute flux calibration. In Figure 2, we have scaled the spectrum to its J 
magnitude assuming a power-law continuum slope of a u = —0.5 (f v ~ v av ). 



3. SPECTRAL FIT AND EMISSION-LINE MEASUREMENTS 

To analyze broad emission lines such as C iv and Mg n we fit and subtract the UV Fe n emission 
that contaminates most of the UV-optical spectral region; thousands of Fe n lines are emitted at 
UV-optical e nergies that blend with othe r lines in the spectrum. We use the UV Fen template 



presented by IVestergaard &: Wilkesl (|200ll ) and adopt the fitting procedure described in detail in 
§ 4.2 by these authors. Briefly, the fitting is done as an iteration over a power-law fit to the 
continuum emission and a model fit to the Fell emission using a scaled and broadened version of 
the Fen template. The Fen template is broadened by convolving the template with Gaussians 
with a range of sigmas. In the first iteration a power-law is fitted to regions with very little or 
no contribution from line emission (e.g. 1275-1295 A; 1325-1335 A; 1425-1480 A; 2180-2220 A; 
3000-3040 A). Upon subtraction of this primary continuum fit multiple broadened copies of the 
Fen template are scaled to regions which predominantly contain Fe emission. The broadened, 
scaled Fen template that provides the best overall match to the observed Fell emission in the 
observed spectrum is then selected and subtracted. Some adjustments of the Fen fitting regions 
are sometimes necessary to optimize the fit, as judged by a minimization of the residual flux upon 
subtraction of the Fe n model. After subtraction of the best matching Fe n model from the original 
data (i.e., with the continuum not subtracted) another power-law fit to the continuum emission is 
performed; the lack of (most of the) contaminating Fe n enables the fitting to be done over a broader 
wavelength range, often resulting in a better and more realistic continuum fit. After subtraction 
of this continuum model, the Fen fitting is repeated. The iteration of the continuum and Fen 
fitting is repeated until these individual fits do not change and the subtraction of the Fen model 



3 The naming convention for SDSS sources is SDSS JHHMMSS.SSiDDMMSS.S, and the positions are expressed 
in J2000.0 coordinates as given in Table 1. We use SDSS JHHMMiDDMM for brevity. 
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leaves a smooth, featureless spectrum (barring the strong, broad non-Fe emission lines) with a 
realistic power-law slope. Convergence is often obtained in the second or third iteration. For SDSS 
J1623+3112, the wavelength range of the spectrum is too short to allow a reasonable fit to both 
the Fen emission and a power-law continuum, so we assume a slope of a u = —0.5 during the model 
fitting. This gives a lower limit of the Fe II em ission. Balmer continuum could be non-negligible at 
A > 2000 A (e.g. Dietrich et alJ 12002 . l2003bl ). however, the wavelength coverage of our spectra is 
not long enough to provide a realistic fit to it, so we do not consider it in our fits. The power-law 
continuum fits to each of the quasars in our sample are listed in column 3 of Table 2. We note that 
the accuracy of the derived line width depends on the quality of the spectra. For widths broader 
than 5000-6000 km/s our ability to accurately determine the intrinsic Fen line width decreases. 

We use Mgn emission lines to determine the redshifts of the quasars. Accurate redshifts are 
important for various studies such as molecular line searches. The redshifts given in the quasar 
discovery papers were mostly determined from the Lya emission lines, where strong Lya absorption 
systems usually cause large uncertainties in redshift measurements. While most of high-ionization, 
broad emission lines such as Civ can be significantly blueshifted with respect to quasar system- 
atic redshifts, the Mg II emission line has a small blueshift and thus provides a relia ble redshift 
measurement (e.g. iMcIntosh et al.lll999l ; IVanden Berk et alj|200ll ; iRichards et al.ll2002l ). TheMgn 
emission line at z ~ 6 is close to a strong telluric absorption band at 1.9 /im, as seen in Figure 
1. To avoid the effect of the telluric absorption, we measure the line center of each Mgn emission 
line using a Gaussian profile to fit the top ~ 50% of the line. The line center is also estimated 
by calculating the 3 x median-2 x mean mode of the top ~ 50% of the line. The two methods give 
similar results because the top ~ 50% of the line is almost symmetric. The measured r edshifts are 
given in Column 2 of Table 2. They are in good agreement with the values derived in iKurk et al 



(|2007l ). but differ by 0.02-0.03 from the values in the discovery papers. We did not detect the 
Mgn emission line for SDSS J0836+0054 and SDSS J1044— 0125, which were observed with high 
humidity. In Tabl e 2, the red s hift o f SDSS J0836+0054, measured also using the Mgn emission 
line, is taken from lKurk et al.l (|2007l ). and the redshift of SDSS J1044— 0125 is measured from the 
C in] A1909 emission line using the same method. 

We fit individual emission lines after the subtraction of the Fe n emission and power-law con- 
tinuum. We use one Gaussian profile to fit most isolated lines and use multiple Gaussian profiles to 
fit blended lines. We allow both the line center and the line width to vary unless stated otherwise. 
Broad emission lines in quasars usually show asymmetric velocity profiles, so in the cases in which 
a single Gaussian profile does not provide a reasonable fit, we use double Gaussian profiles. There 
is no particular physical meaning to the two Gaussian components. For each quasar, the following 
emission lines (if detected) are fitted. (1) Lya + Nv + Sin A1262 (hereafter Sin). We use four 
Gaussian profiles to fit the three lines, with the first two profiles representing the broad and narrow 
components of Lya and the last two profiles r epresenting Nv and Sill, r espectively. During the fit, 
the line centers are fixed at the values given by lVanden Berk et alJ (|200ll ). Since the blue side of the 
Lya line is strongly absorbed, we only fit the red side of the Lya line. (2) O I A1304, C n A1335, and 



-7- 



Si IV A1397 (hereafter O I, C II, and Si iv). We use one Gaussian profile to fit each of these lines. (3) 
Civ. The strong Civ line is fitted using one or two Gaussian components. (4) Hen + Oiii] Af663 
(hereafter Oiii]). These two lines are very weak, but Hen is important for determining chemical 
abundances. They were tentatively detected in three quasars and fitted simultaneously using two 
Gaussians. (5) Aim A1857 + Cm] A1909 (hereafter Aim and Cm]). The weak Aim line is often 
detected on the blue wing of the strong Cm] line. We use two Gaussian profiles to fit these two 
lines simultaneously if Aim is detected, otherwise the Cm] line is fitted using a single Gaussian 
profile. We ignore the weak Sim] A1892 line on top of Cm], (6) Mgll. The Mgn emission line 
can be well fitted by double Gaussians (e.g. iBarth et al.ll2003l ). with a weak component located 
on the blue side of a much stronger component. If the blue side of Mgn is affected by telluric 
absorptions, we fit the line using a single Gaussian. For the emission lines that are severely affected 
by broad-absorption lines (BALs), we fix their central wavelengths and fit their red half sides using 
a single Gaussian profile. 

Figures 1 and 2 show the fitting results for the six quasars. The blue solid lines are the 
best fits to the emission lines and the blue dashed lines are the best model fits to the power-law 
continuum and Fe n emission. The red lines are the sums of all the compo nents. For comparison, 
we show the low-redshift composite spectrum of IVanden Berk et al.l (|200ll ) in gray. Table 3 shows 
the rest-frame EW and FWHM of emission lines derived from the model fitting results. For the 
lines which were fitted using double Gaussian profiles, their EW and FWHM are calculated from 
the combination of the two components. The errors given in the table are the formal uncertainties 
obtained from our fitting process, so the real errors could be larger in the cases in which the lines are 
heavily blended, or we fit a profile of the wrong shape. The He n emission is tentatively detected 
in three quasars at a significance level of ~ 2a. We note from Table 2 that three quasars have 
continuum slop es consistent with the slo pe a u = —0.44 of the composite spectrum of low-redshift 
SDSS quasars (|Vanden Berk et alJhoOl h. w hile SPSS J1030+0524 and SDSS J1306+0356 show 
significantly bluer continua with a u ~ 0.5. iPentericci et al.l (|2003l ) obtained an average slope of 
a u = —0.57 with a scatter of 0.33 from 45 intermediate-redshift SDSS quasars at 3.6 < z < 5.0. 
This indicates that the continuum slopes in the two unusual quasars deviate from the average slopes 
at a ~ 3a significance level. Table 2 also shows the rest-frame op tical continuum s lopes derived 
from the Spitzer IRAC broad-band photometry of the same objects (| Jiang et al.ll2006l ). The optical 
continua ar e generally redder than th e UV continua measured from the Gemini spectra, as also 



reported by lVanden Berk et al 



(2001 



4. CHEMICAL ABUNDANCES AND CENTRAL BH MASSES 



4.1. Chemical Abundances 



Emission-line ratios can be used to measure gas metallicity in the quasar BLR and track 
the chemical evolution with redshift. When weak emission lines are difficult to detect in high- 
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redshift quasars, strong broad emission lines such as N v, C IV, and Mg II become a powerful tool to 
probe metallicity at high redshift (e.g. lHamann Ferlandlll992l . Il993l ). Different elements form on 
different timescales. Elements such as C, O, and Mg are formed in the explosions of massive stars, 
and their enrichment is rapid; while the second generation element N is produced from previously 
generated C and O, and its enrichment is relatively slow. Helium is a primordial element and its 
abundance changes little with cosmic time. Chemical abundances measured from emission lines 
and intrinsic absorption lines sho w solar or supersolar metallicity in quasar enviro nments at low 



redsh ift (e.g. lHamann et al.ll2007l ). and can be as high as 15 times solar metallicity ([Baldwin et al 
2003). Furthermore , studies haye sho wn that there is little chemical evolution up to the highest 
known redshift (e.g. iBarth et al.l l2003). 



We calculate emissiondine fluxes from the model fitting results described in § 3. The fluxes 
are normaliz e d to the Civ fluxes and shown in Table 4 and Figure 3 (filled circles with errors). 
Nagao et al.l (|2006l ) used a sample of over 5000 quasars from the SDSS Data Release Two to 
study quasar BLRs. They made quasar composite spectra in the ranges of 2.0 < z < 4.5 and 
—29.5 < Mb < —24.5, and measured emissiondine ratios and metallicities in the composite spectra 
for each redshift and luminosity bin. They found that most of the line ratios (with respect to the 
C iv line) that they investigated do not show strong evolution with redshift. For the purpose of 
comparison, we show in Figu re 3 the line ratios of quasars (filled squares) in the luminosity range 
of —27.5 < Mb < —26.5 from lNagao et al.l (120061). A lthough there are correlations between quasar 
luminosities and most line ratios (jNagao et al.l 120 06), most of our quasars have luminosities in a 
small range of —27.5 < Mb < —26.5. Figure 3 shows that most of the flux ratios do not exhibit 
strong ev olution up to z ~ 6 . There is a trend for higher Cii/Civ and Oi/CiV ratios at higher 



redshift. iNagao et al.l (120061 ) already noted that the Oi/Civ flux ratio is marginally correlated 
with redshift. This trend could also be caused by the intrinsic dispersion of the flux ratios. 

Photoionization models have shown that a series of emissiondine ratios can be used to estimate 



gas metallicity in the BLRs of quasars. lHamann et al.l (|2002l ) studied the relative N abundance as 



a me tallicity indicator based on the locally optimally emitting cloud (LOC) model (IBaldwin et al 



19951 ) , and calculated theoretical emissiondine ratios as a function of metalli city Z/Zq. We e s timat e 
the metallicity from the Nv/Civ and Nv/Heii ratios using Fig ure 5 o f Ham ann et al.1 (J2002j). 
There is a 30% differen ce between the sola r abundances used by Hamann et al.l ( 20021') and the 



latest solar abundances ([Baldwin et al.l l2003). so we scale the abundan ces of lHamann et al.1 ((2002) 
to match the latest values u sing the method by IBaldwin et al.l (|2003l ) . The results are shown as 
the open circles in Figure 4. INagao et al.l (|2006l ) also used the LO C model, but i nclud ed more UV 
emission lines in their photoionization model. We use Figure 29 of lNagao et al.l (|2006l ) to measure 
the metallicity from various line ratios and show the measurements as the filled circles in Figure 4. 
The two models give similar metallicities from the Nv/Civ and Nv/Heii rati os. The meta ll icities 



in these high-redshift quasars are supersolar, with a typical value of ~ 4 Z Q . INagao et al.1 (J2006) 
computed chemical abundances as a function of redshift for the SDSS composite spectra using 
their photoionization model. The filled squares in Figure 4 show the abundances averaged in the 
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luminosity range of —28.5 < Mb < —25.5. The fill e d trian gles represent the chemical abundances 
of eleven 3.9 < z < 5.0 quasars from iDietrich et al.l (|2003al ). Figure 4 shows that the metallicity is 
consistent with no strong evolution within the errors up to z ~ 6. 

The Fe n /Mg n ratio is important in understanding chemical evolution at high re dshift. We use 
the Fe II emission line co mplex at Aq = 2000 ~ 3000 A as the Fe n flux indicator (e.g. IDietrich et al 



2002 : iBarth et al.1 12003 ). The flux of the Fen complex is integrated from 2200 to 3090 A over 
the best-fitting Fen template. The derived Feii/Mgii ratios (filled ci rcles) with a typic al value of 
4.9 ±1.4 are shown in Table 4 and Figure 5. Using a similar method iKurk et al.l (|2007l ) obtained 
a lower ratio of 2.7 ± 0.8 from their quasar sample. The discrepancy is partly due to the fact 
that they corrected for Balmer continuum in their spectral fits. For comparison, we also show the 



Feii/Mgii ratios of z ■ ~ 6 quasars from previou s measurements (jBarth et al 



2003 



Freudling et al 



2003 



Iwamuro et al.l 12004 ) as open circles in Figure 5. 



2003; 



Maiolino et al 



Iwamuro et al 



measured the Fell/Mgn ratios from the spectra of quasars at < z < 5.3, and computed the 
median values of the ratios for a range of redshifts. These median ratios are show n as the filled 
squares in Figure 5. To study the evolution of Feii/Mgii, IDietrich et al.1 (|2003bl ) made quasar 
composite spectra from z = to 5. The Feii/Mgii ratios measured from the composite spectra are 
shown as the filled triangles in Figure 5. We note that the Fen flux in some studies was integrate d 



over the wavelength range of 2150 < Aq < 3300 A (e.g. iFreudling et al.ll2003l ; llwamuro et al.ll2004l ) 



slightly different from the range that we used. We also note that some of these studies used different 
Fen templates to measure Fen fluxes. These issues usually affect the Feii/Mgii measurements by 
less than a factor of two and cause a relatively large scatter in the measurements at z ~ 6. Figure 5 
shows that our results are consistent with those derived from both low-redshift samples and other 
z ~ 6 samples within errors. Becau se the details of how the Fen bump is formed are not well 
understood (e.g. Baldwin et al.l 120041 ) . so we do not derive actual abundances from the Feii/Mgii 
ratios. 

We have shown that the metallicity in the BLRs of high-redshift quasars is supersolar, and the 
lack of strong evolution in metallicity continues to z ~ 6. The high metallicity at z ~ 6 indicates 
that vigorous star formation and element enrichment have occurred in quasar host galaxies in the 
first Gyr of cosmic time. Millimeter and submillimeter observations revealed that luminous z ~ 6 
quasars are ext remely far-infrared (FIR) luminous (~ 10 13 Lm) and have a large amount of cool d ust 
(10 8 -10 9 M ) dCarilli et alJboOll : iBertoldi et alJl2003al : IPriddev et alJl200.4 Ifiobson et aOhood ). If 
the FIR emission is mainly powered by star formation in host galaxies, the star formation rates are 
estima ted to be ~ .1000 Mm yr -1 . Stron g starbursts can be induced by merging gas-rich galaxies 



(e.g. iHopkins et al.ll2006l ; iLi et al.1 120071 ) . CO observations have already revea led the presence of 
~ 2 x 10 10 Mm of molecular gas in the hig hest redshift quasar SDSS J1148±5251 (jWalter et alJl2003l : 
Bertoldi et al.l l2003bl ; IWalter et al.l 12004 ) . This amount of gas can sustain the star formation rate 
of ~ 3000 Mq yr _1 , inferred from the FIR lu minosity, for 10 7 -10 8 years. If the typical duty cycle 



time of luminous quasars is 10 7 -10 8 years (e.g. iKauffmann Haehneltl 



2000 



Wvithe k Loeb 



2003 



Shen et al.1 120071 ) . it would have enough time to form a large amount of first generation elements 
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from massive stars. 



We have also shown that the Feii/Mgn ratios of the z ~ 6 quasars are consistent with low- 
redshift measurements. Most of the Fe in the solar neighborhood is produced in intermediate-mass 
stars, and their enrichment is delayed by ~ 1 Gyr compared to Mg. Mg + and Fe + have similar 



ionizat ion potentials and thus the Fen/Mgn ratio reflects the Fe/Mg abundance (IDietrich et al 



2003bl ). Therefore, if Fe was generated in the same way as it is in our neighbors, the Fen/Mgn 
ratio is expected to be a strong function of cosmic time at high redshift. However, we do not 
see any evolution in this ratio even at z ~ 6. It has been shown that the timescale of Fe enrich- 
ment from SNe la strongly depends on star formati on histories of host galax i es, and can be much 



short er in high-redshift qu a sar en vironments (e.g. iFriaca fe Terlevich 



1998; Matteucci fc Recchi 



2001 



). iMatteucci fe Recchil (120011 ) found that the timescale for the maximum SN la rate (Fe en- 
richment) can be as short as ~ 0.3 Gyr in an ellip tical galaxy with a high star formation rate but 
a short star formation history. On the other hand. IVenkatesan et al.1 (120041 ) pointed out that stars 
with a present-day initial mass function are sufficient to produce the observed Fe n /Mg n ratios in 
z ~ 6 quasars, and SNe la are not necessarily the main contributor. Fe could also be generated in 
Population III stars, which are suggested to be first generation stars with masses between 100 and 
1000 Mq. The se very massive stars p roduce a large amount of Fe within a few Myr due to their 



short lifetimes (jHeger fe Wooslevl 120021 ) . 



4.2. Central BH Masses 



The central BH masses of high-redshift quasars are important in understanding the growth 
of BHs and quasar accretion rates in the early universe. For high-redshift quasars where direct 
mass measurements are difficult, BH masses can be estimated using mass scaling relations based on 
broad emission line widths and continuum luminosities. S trong emission lines , including H/3, Civ, 



and Mg II, have been used to determine BH masses (e.g. IWandel et al.1 Il999i : IDietrich fe Hamann 



20041 : iMcLure fe Dunlonlbooi Ivestergaard fe Petersonlbood l . Based on empirical scaling relations, 
a few luminous z ~ 6 quasars have been fou nd to harbor central BHs with masses of 10 9 ~ 10 10 



Mq (e.g. lBarth et al.1 12003 : IVestergaardl 120041 1. 



We estimate the central BH masses for the six quasars using the following BH mass scaling 
relations, 

/FWHM(Civ)\ 2 / AL A (1350A) V 



Mbh(Civ) =4.57 



V km s 



0.53 



V10 44 ergs s- 1 J 



M 







and 



\ /FWHM(Mgn)\ 2 / AL A (3000A) \°' 62 _ 
M BH (Mgii) = 3.2f km { s J j ) ( „„;, , ) M 



\10 44 ergs s- 1 J 



(1) 



(2) 



derived by IVestergaard fe Peterson! (|2006l ) and lMcLure fe Dunlopl (|2004j), respectively. L A (1350A) 
and L A (3000A) are luminosities at rest-frame 1350 and 3000 A. The results are listed in Columns 
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4 and 5 of Table 5. They are in agreement with the masses derived by iKurk et al.l (|2007l ). For 
SDSS J1030+0524 and SDSS J1306+0356, the Civ-based BH masses are greater than the Mgn- 
based BH masses by a factor of ~3, but still within the intrinsic scatter (a factor of 3-4) in these 
scaling relations. The unusual blue continua in the two quasars could cause this discrepancy. The 
Mgn-based BH masses are also estimated using the latest scaling relation based on FWHM(Mgn) 
and La(1350A) by Vestergaard et al. (in preparation). The results are shown in Column 6, and 
they are consistent with the ClV-based BH masses. Similar to other luminou s quasars at z ~ 6, 



our quasars have central BH masses of 10 9 -10 10 M (e.g. Earth et al 
Jiang et alJhood ). 



2003 



Vestergaard 



2004 



We calculate the Eddington luminosity ratios of the quasars using their bolometric luminosities 
and BH masses. We use the Mgn-based BH mass for SDSS J1623+3112 and the Civ-based BH 
masses for the others. The bolom etric luminosities were calculated from multiband observations of 
the quasars by I Jiang et al. The results are given in Column 3 of Table 5. These luminous 

z ~ 6 quasars have E ddington ratios of order unity, comparable to quasars with similar l uminosities 



at lower redshift (e.g. iMcLure fc Dunlopll2004l : IVestergaard 



2004 



Kollmeier et al.ll200fih . 



It is remarkable that billion-solar-mass BHs can form less than one Gyr after the Big Bang. 
With the reasonable assumption that BH accretion is at the Eddington rate, the BH mass at time 
t is Mt = Moe*^ 1 ", where e is the radiative efficiency, r = 4.5 x 10 8 years, and Mq is the initial BH 
mass or the seed BH mass. Seed BHs can be produced fro m the collapse of Population HI stars or 



gas c l ouds, and their masses are roughly 10 2 -10 4 M (e.g. iMadau &: Rees 



2001 



Volonteri fc Ree; 



2006 



Lodato Natarajanl 120071 ). Consider the case in which a massive BH at z = 6 formed from 
a seed BH with Mq = 10 3 M at z = 20. The e-folding time er for the BH growth is roughly 
4.5 x 10 7 years if e ~ 0.1. The BH grows from z = 20 to 6 by 15 e-foldings, or a factor of 3.3 x 10 6 , 
which results in a massive BH with Mt = 3.3 x 10 9 M at z = 6, comparable to the observed 
BH masses in our sample. If a quasar is shining at half of the Eddington limit, its BH grows 
from z = 20 to 6 by only 7.5 e-foldings, or a factor of ~ 2000, making it very difficult to form 
billion-solar-mass BHs in this scenario. In addition, if Eddington-limited accretion is via standard 
thin disk s, BHs are likely to be spun up and the radiat ive efficiency and Eddington timescale will 
increase (jVolonteri &; Reesll2006l : iRees fc Volonterill2006l ). In this case it would take much longer to 
form massive BHs. So super-Eddington accretion or lower radiative efficiency is probably required 
to form BHs with M t = 10 9 -10 10 M by z = 6. 



4.3. Notes on Individual Objects 



SDSS J0836+0054 (z = 5.810). SDSS J0836+0054 was discovered bv iFan et al.1 (|200lh . 
Its redshift estimated from the Lya emi ssion line is 5.82. comparable to the redshift 5.810 ± 0.003 
determined from the M g II emission l ine (IKurk et al.ll2007l ) . Due to high humidity, we did not detect 
the Mgn emission line. IStern et al.1 (|2003l ) found a very red continuum slope of a v = —1.55 in their 
NIR spectrum of this quasar. However, using the GNIRS spectra with a longer wavelength coverage 
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and a higher spectral quality we obtaine d a slope of a 



slope of the SDSS composite spectrum (jVanden Berk et al.l 12001 



0.62, which is close to a u = —0.44, the 
SDSS J0836+0054 is the most 



luminous quasar known at z > 5.7. The central BH mass estimated from the Civ emission line is 



Becker et al. 


1995) 


Petric et al. 


(2003) 



1995). The flux at 1.4 GHz measured by FIRST is 1.11±0.15 mJy and measured by 



(|2003l ) is 1.75±0.04 mJy. 



SDSS J1030+0524 (z = 6.309). SDSS J1030+052 4 is the second highest redshift quasar 
known to date. The redshift estimated from Lya is 6.28 (|Fan et al.ll200ll ). We obtained a more 
accurate and s l ightly higher redshift 6.309 ± 0.009 from Mgn, consistent with the value 6.311 of 
Iwamuro et all (|2004l ) . The emission- line ratios and m etallicity in this object have been studied by 
a few groups. In the discovery paper iFan et all (|200ll ) estimated lower limits for t he N v/C IV and 



N v/H e II ratios: N v/C iv > 0.4 and Nv/He n > 3.0. Using VLT NIR observations IPentericci et al 
(|2002l ) found Nv/Civ = 0.35 and Nv/Heii > 4.3. From our measurements we obtained Nv/Civ 
= 0.21 ±0.02 and Nv/He II = 2.85±1.70, comparable to the previous studies. The F ell/Mgll ratios 



of this object measured i n prev ious papers are quite differe nt. The values reported by lMaiolino et al 



(12003! ) . iFreudling et all (|2003l ). and llwamuro et all (|2004! ) are 8.65 ± 2.47, 2.1 ± 1.1 and 0.99+/, 



.+1.86 



respec tively. We obtained Feii/ Mgil = 5.46 ± . 90, co nsist ent with the re s ult of iMaiolino et al 



(|2003l ). Note that the Fell flux in IFreudling et al.1 (|2003l ) and llwamuro et all (|2004l ) was integrated 
over th e wavelength range 21 50 < Ao < 3300 A, slightly different from the range 2200 < Ao < 3090 
A that IMaiolino et al.1 ( 20031 ) and we used. In addition, our NIR spectra have much higher SNRs 
than those of previous studies. 



SPSS J104 4-0125 (z = 5.778). SDSS J1044-0125 was the first q 



(|Fan et al.1 12000| ) . The redshift estimat ed from Lya is 5.80 (|Fan et al 



from Civ is 5.74 (Goodrich et al 



2001 



uasar discovered at z > 5.7 



2000) and that estimated 



H owever, these measur e ments could be severel y affected 



by absorption since this is a BAL quasar ((Goodrich et al 



2001 



Djorgovski et al.ll200ll ). Due to 



high humidity we did not detect its Mg n emission line. The redshift measured from C in] (Figure 
1 shows that this line is not affe c ted b y absorption) is 5 . 778 ± 0.005, consistent with the value 
5.78 obtained by IFreudling et al.1 (|2003l ). IFreudling et al.1 ([2003) also reported a Fen/Mgn ratio 
of 5.0 ± 2.1 for this quasar. SDSS J 1044— 0125 has the highest BH mass of the ob jects in our 



samp le. Submillimeter observations (jPriddey et al.1 l2003l ) and Spitzer observations (IJiang et al 



2006) revealed a large amount of cool and hot dust in this object. SDSS J1044— 0125 is the only 
BAL quasar in our s ample. The fracti on of BAL quasars in the sample is 16.7%, similar to the 
low-redshift fraction ([Trump et a.l 120061 ) . 



SDSS J1306+0356 (z = 6.016). The d iscovery paper for SDSS J1306+0356 indicated a 
redshift z = 5.99 based on Lya (IFan et al.ll200ll ). We obtained a more accurate and slightly higher 



redshift 6.016 ± 0.005. IPentericci et al.l (120021 ) reported a Nv/Civ ratio of 0.67 for this quasar. 
The ratio we measured is 0.55 ± 0.05. We al so obtained a Fell/M gii ratio of 5.77 ± 0.72, slightly 
smaller than the value 9.03 ± 2.26 derived by lMaiolino et al.1 (120031 ) based on low SNR spectra. We 



detected four strong Mgn AA2796,2803 absorption systems in the spectra of this quasar, including 
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the two highest redshift Mgn absorbers known (see § 5). 



SDSS J1411+1217 (z = 5.927). SDSS J1411+1217 has the narrowest emission lines in our 
sample. Its Lya and N v emission lines are well separated. Therefore the redshift z = 5.93 measured 
from Nv and Ol in the discovery paper (IFan et al.l 12004 ) is accurate, and is consistent with the 
redshift 5.927 ± 0.004 measured from Mgil in the Gemini spectrum. SDSS J1411+1217 has the 
lowest BH mass in the sample. It also has a very weak Spitz er 24/j.m (Ao ~ 3.5 fim) flux, and its 
SED does not show any hot dust emission (| Jiang et al.l 120061 ) . Another quasar SDSS J0005— 0006 
(not in this sample) shows similar IR properties and does not have hot dust emission either. This 
may reflect dust evolution at high redshift, as such dustless quasars are unknown at lower redshift. 

SDSS J1623+3112 (z = 6.247). The redshift of SDSS J1623+3112 measured from Lya in 
the discovery paper is 6.22. We obtained a higher redshift of 6.247 ± 0.005 from Mgil. This makes 
it the third highest redshift quasar known to date. 



5. Mg II ABSORPTION SYSTEMS AT z < 6 

We search for Mgil AA2796,2803 intervening systems in the spectra of five z ~ 6 quasars. SDSS 
J1623+3112 is excluded in this analysis because it has a short wavelength coverage. The SNRs of 
the spectra for different quasars are different, as are the dispersions from one order to another. The 
SNRs at A < 1.2 fim are usually higher than those at longer wavelengths, especially in the range of 
1.4 < A < 1.8 /um (roughly H band). However, a Mgil A2796 absorption line with rest equivalent 
width Wq 2796 greater than 1.0 A at z < 3 or a Mgil A2796 absorption line with Wq 2796 greater 
than 1.5 A at z > 3 can be detected at a ~ 7a significance level in the lowest SNR parts of the 



spect ra shown in Figure 1. Moreover, previous studies have used similar limits (e.g. iNestor et al 



20051 ). Our sample is small, and does not provide good statistics by itself, so the comparison with 
other samples is important. Thus to achieve unbiased statistics, we search for Mgil absorbers with 
W A2796 > 1.0 A at z < 3 and Mgil absorbers with W A2796 > 1.5 A at z > 3. 

We search for pairs of absorption lines with separation matching the Mgil AA2796,2803 doublet. 
Then we search for absorption lines of other species associated with the identified Mg n absorbers 
by matching various strong quasar absorption lines known in the literature to the spectra. Five 
Mgil absorption systems are identified, with one in the spectrum of SDSS J0836+0054 and the 
other four in the spectra of SDSS J1306+0356. We normalize the absorption lines using the local 
continua measured by fitting low-order spline curves to the spectra. Figure 6 shows the normalized 
spectra of the five absorption systems, with the identified lines marked. These lines are usually 
strong and commonly detected in quasar absorption systems, such as Civ AA1548,1550 (marked as 
Civ-D in the figure), Fen A1608, Alii A1670, Aim AA1854,1862 (marked as Alm-D), Fen A2600, 
Mgil AA2796,2803 (marked as Mgn-D), and Mgi A2852. The four absorption systems found in the 
spectra of SDSS J1306+0356 are denoted as (a), (b), (c), and (d). 

We measure the redshifts (z a b s ) and rest-frame equivalent widths (Wo) of the absorption lines 
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by fitting a Gaussian to each line. Unresolved or very weak absorption lines are not fitted. The 
results are given in Table 6. The Mg n doublet ratios are close to one, indicating that these strong 
Mgn absorbers are saturated. The Mgn absorber identified in the spectra of SDSS J0836+0054 
is at z = 3.742. A pair of moderate strong Aim AA1854,1862 lines in the same absorber are also 
detected. C iv AA1548,1550, which is usually associated with Mgn, is beyond the spectral coverage. 
We do not detect any strong Fe n absorption in this system. There are four Mg n absorption systems 
identified in the spectra of SDSS J1306+0356. Systems (a) and (b) are at relatively low redshifts. 
System (a) at z = 2.202 has a strong Mgi A2852 absorption line. System (b) at z = 2.532 has 
a strong Fell A2600 absorption line. Civ AA1548,1550 in both systems is beyond the spectral 
coverage. We discovered two strong Mgn AA2796,2803 absorbers at z = 4.8668 and 4.8823. They 
are the most distant Mg n absorbers known to date. The Mg n A2803 line in system (c) and the Mg n 
A2796 line in system (d) almost exactly overlap with each other, so we do not calculate Wq for the 
two lines. Note that this is not the so-called line-locking effect , which mainly occurs in BAL quasars 
or other associated absorp tion systems with z a b s ~ z em (e.g. iFoltz et al.lll987l ; IVilkoviskii &; Irwin 
200ll : iFechner et al.l 12004 ) . since the two absorbers are far from SDSS J1306+0356. The proper 
distance between the two absorbers is about 1.4 Mpc, so they are not likely associated with the 
same galaxy, but could still belong to the same large-scale structure. Both absorption systems 
have strong Alii A1670 and Fen A2600 absorption and weak Fen A1608 and Aim AA1854,1862 
absorptions. The Civ AA1548,1550 absorption is strong in system (c) but very weak in system 
(d). It has been found that Civ AA1548,1550 is alm ost always detected in strong Mgn absorption 
systems at low redshift (e.g. ISteidel Sargentlll992l). The wea kn ess of Civ in s ystem (d) indicates 



that it is a rare, C iv-deficient system ([Churchill et al.l 12000 ) . lElston et al.l (|1996l ) discovered a 
strong Mgn absorber at z = 4.38 in the spectrum of a z = 4.68 quasar. They found that the Civ 
absorption in this system is also very weak. Thus among the three strong Mg n absorbers at z > 4, 
two have very weak associated C iv, probably indicating t hat z > 4 absorpt ion systems are in a 



lower ionization state than their low-redshift counterparts (jElston et al.lll996l ). 



We calculate the comoving line-of-sight number densities (dN/dz) of the Mgn absorbers. The 
sample of the five Mgn absorbers is divided into two subsamples, one with W A2796 > 1 A at z < 3 
and the other one with W 7 ^ 2796 > 1.5 A at z > 3. The redshift path covered by a spectrum is, 



AZ(iy A2796 ) 



g(W X2796 ,z)dz, 



(3) 



where g(W^ 



A2796 



1 if a Mgli absorber at redshift z and of the given EW would be detected 



in the spectrum in the redshift range Az and ^(Wq 42796 , z) = otherwise. The total redshift paths 
covered by the two subsamples are 3.4 and 10.4, respectively. Then the number density along the 
redshift path for each subsample is estimated by, 

N a 



dN/dz 



£ AZ(^ 2796 ) 4 



(4) 



where N a is the number of Mg n absorbers and the sum is over all spectra in the subsample. We 
show our results in Figure 7. The filled diamond represents the density for Wq 12796 > 1.0 A at z < 3 
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and the filled circle represents the density for W^ 2796 > 1.5 A at z > 3. The two data points are 
positioned at the median redshifts z = 2.6 and 4.3, respectively. The statistical uncertainties are 
Poisson. 



The number density of Mg n absorbers at low redshift has been well studied. ISteidel Sargent 



(|1992l ) collected the spectra of 103 quasars and found that the density of Mgn absorbers with 

pj/A2796 > q 3 ^ doeg not showg 

c osmic evoluti o n, bu t the density of strongest Mgn absorbers sig- 
nificantly increases with redshift. iNestor et al.l (|2005l ) reached similar results by st udying Mg n ab- 



sorpt ions in the spectra of thousands of quasars from the SDSS Early Data Release (jStoughton et al 



2002 ). They divided the sample into subsamples with W A2796 > 0.3, 0.6, 1.0, 1.5, 2.0, 2.5, and 3.0 
A, respectively, and concluded that dN/dz of the Wq 2796 > 0.3, 0.6, 1.0, and 1.5 A subsamples are 
consistent with no evolution within errors, while the densities of the other three subsamples show 
strong evolution with redshift. For co mparison, we show in Figure 7 the W i V 2796 > 1.0 and 1.5 A 
samples (open diamonds and circles) of lNestor et al.1 (|2005l ). iProchter et al.1 (120061 ) studied strong 
Mgll absorption system s using an even larger quasar sample from the SDSS Data Release Three 



( Abazajian et al.ll2004l ). They found that there is no evolution in the densities of W^ 2796 > 1.0 A 
absorbers at 0.8 < z < 2.0. The crosses in Figure 7 show the densities for their Wn 279e > 1.0 A 



sample (they did not use a sample of W A2796 > 1.5 A). In Fi gure 7 we also show the non-evolution 
curves (dashed lines) for the cosmology of (Q m ,Q\, h) = (0.3,0.7,0.7) with zero curvature, 



dN/dz = N (l + z) 2 [£l m {l + zf + 0/ 



-1/2 



(5) 



w here No is a normal ization constant. The curves have been scaled to minimize the \ 2 to the data 



of lNestor et al.l ([20051 ) . The density of Mgn absorbers with W A2796 > 1.5 A that we have measured 



is consistent with no cosmic evolution up to z > 4. However, due to the large uncertainty, we 
cannot rule out the density evolution of Mgn absorbers. 



6. SUMMARY 

We have obtained high-quality NIR spectra of six luminous z ~ 6 SDSS quasars from our 
Gemini observations. Five of them were observed using Gemini-South/GNIRS in cross dispersion 
mode, which provides a simultaneous coverage of 0.9-2.5 (im. The sixth was observed in K band 
with Gemini-North/NIRI. These spectra have higher SNRs and better wavelength coverage than 
those in previous studies. We use the NIR spectra combined with optical data to study chemical 
abundances and BH masses in the z ~ 6 quasars and strong intergalactic Mgn absorption at 
2.2 < z < 6. 

The spectra are fitted using a combination of a power-law continuum, an Fe n emission tem- 
plate, and a series of emission lines. We calculate the fluxes for detected emission lines based on 
the Gaussian model fitting results. We find that the line flux ratios (normalized to the Civ fluxes) 
that we have investigated do not show strong evolution with redshift. We calculate gas metallicity 
from emission-line ratios using the results of photoionization models given by previous studies and 
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find that the metallicity in the BLRs of these high-redshift quasars is supersolar with a typical 
value of ~ 4 Z . The comparison with low-redshift observations shows no show strong evolution in 
metallicity up to z ~ 6. The Feii/Mgii ratio is also measured for each quasar. We find a typical 
Fell/Mgll ratio of 4.9 ± 1.4, which is consistent with low-redshift samples. All these measurement 
indicate the existence of vigorous star formation and element enrichment in host galaxies in the 
first Gyr after the Big Bang. We have estimated central BH masses from the C iv and Mg n emis- 
sion lines using empirical mass scaling relations. As found in other luminous z ~ 6 quasars, these 
quasars have BH masses of 10 9 -10 10 M© and Eddington luminosity ratios of order unity. 

We have searched for strong Mgn AA2796,2803 intervening systems in the spectra of five 
quasars. Two Mgli absorbers with rest equivalent width Wq 279 ® > 1 A at 2.2 < z < 3 and three 
absorbers with Wq 279 ® > 1.5 A at z > 3 are identified in the spectra of two quasars. The two 
most distant absorbers are at z = 4.8668 and 4.8823, respectively. These are the highest redshift 
Mgli absorbers known to date. We calculate the comoving line-of-sight number densities for the 
five identified Mgli absorbers. By comparing with low-redshift studies we find that the densities 
(dN/dz) of Mgli absorbers with W A2796 > 1.5 A are consistent with no cosmic evolution up to 
z > 4. We note that our sample is small. A larger sample is needed to provide a good constraint 
on the density evolution of Mgli absorbers. 
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Fig. 1. — Optical and NIR spectra of five z ~ 6 quasars in our sample. The spectra at A > 1.0 
/im were obtained with Gemini-South/GNIRS and the spectra at A < 1.0 /jm are taken from the 
quasar discovery papers. The spectra have been smoothed by three pixels. The blue solid lines are 
the best fits to emission lines and the blue dashed lines are the best fits to the power-law continuum 
and Fe n emission. The red lines ar e the sums of all the compo nents. For comparison, we show the 
low-redshift composite spectrum of IVanden Berk et al.1 (|200ll ) in gray. A spectrum of a standard 
star HIP 54027 is also shown in the figure with a black body (with the effective temperature of 
HIP 54027) spectrum (in red) on top of it. The bottom panel shows the sky transparency at an 
airmass of 1.0 and water vapor of 1.6 mm. 
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Fig. 2. — Same as Fig. 1 for SDSS J1623+3112, which was observed in K band using Gemini- 
North/NIRI. 
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Fig. 3. — Emission line fluxes (filled circles with errors) of the z ~ 6 quasars compared to low- 
redshift measurements. The fluxes have been normalized to the C iv fluxes. The filled squares 
represent the flux ratios me asured in the comp osite spectra of quasars in the luminosity range 
—27.5 < Mb < —26.5 from INagao et al.1 ((2006). Most of the flux ratios do not exhibit strong 
evolution up to z ~ 6. 
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Fig. 4. — Chemical abundances of quasars derived from line ratios as a function of redshift. The 
ope n circles represent the metallicities derived from the Nv/Civ and Nv/Heii ratios using Figure 
(200 2]). Th e filled circles represent the metallicities derived from the line ratios 



5 of Hamann et al 



using Figure 29 of lNagao et al.l (|2006l ) . The filled squ ares show the avera ged metallicities of quasars 
for the luminosity range —28.5 < Mb < —25.5 from lNagao et al.l (120061 ) . The filled triangles show 



the metallicities of eleven 3.9 < z < 5.0 quasars from 



Dietrich et al. 



(|2003al ). The metallicities 



estimated from each line ratio show no strong evolution up to z > 6 within the errors. 
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Fig. 5. — The Fen/Mgn ratio as a function of redshift. The filled circles are our results. The open 
circles show the Feil/Mgil ratios of z ~ 6 quasars from previous measurements (jBarth et al.ll2003l : 



Maiolino et al 



2003; 



Freudling et al.1 L2003 ; 



represent the low-redshift results from 



Iwamuro et al 



Iwamuro et al 



2004 ). The filled squares an d triangles 



(|2002l ) and [Dietrich et all (boOSbl ). The lack 
of strong evolution of the Fen/Mgn abundance continues to z ~ 6. 
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Fig. 6. — Normalized spectra of the five absorption systems in the spectra of SDSS J0836+0054 
and SDSS J1306+0356. The positions of the identified lines are marked as vertical lines, and the 
names of the lines are given nearby. Doublets are expressed as '— D'. The four absorption systems 
in the spectra SDSS J1306+0356 are denoted as (a), (b), (c), and (d). 
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Fig. 7. — Comoving number densities of Mgn absorbers. The filled diamond represents the density 
for Wq 2796 > 1 A at z = 2.6 and the filled circle represents the density for Wq 2796 > 1.5 A at 
z = 4.3. The open diamonds and circ les show the number densities for Wq 279G > 1.0 and 1.5 
Nestor et al.1 (120051 ) . The crosses show the densities of the Wq 2796 > 1.0 A 
(|2006l ). The dashed lines represent non-evolution curves for the 



A from Figure 9 of 



absorbers from Prochter et al 



cosmolo gy of (Q m ,Q\, h) = (0.3,0.7,0.7). The curves have been scaled to minimize the x 2 to the 
data of lNestor et al.l (|2005l ). The densities are consistent with no cosmic evolution. 



Table 1. Log of Observations 



Quasar fSDSS"> 


Redshift 


y a t> \ in a0 i 
^AB \ LLLcl E>J 




Date 




In «tvn in pnt" 

lllOll LllllClllj 


f ( 777 1 77 I 


J083643.85+005453.3 


5.82 


18.74±0.05 


17.89±0.05 


2006 Feb 19 




GNIRS 


70 


J103027.10+052455.0 


6.28 


20.05±0.10 


18.87±0.10 


2006 Mar 20, 


26 


GNIRS 


140 


J104433.04-012502.2 


5.80 


19.23±0.07 


18.31±0.10 


2006 Mar 28 




GNIRS 


120 


J130608.26+035626.3 


5.99 


19.47±0.05 


18.77±0.10 


2006 Mar 26, 


27 


GNIRS 


140 


J141111. 29+121737.4 


5.93 


19.65±0.08 


18.95±0.05 


2006 Mar 28 




GNIRS 


120 


J162331.81+311200.5 


6.22 


20.09±0.10 


19.15±0.10 


2004 Aug 16, 


Sep 14 


NIRI (K) 


180 










2005 Aug 22, 


24 







Note. — Reds h ifts, zar (AB m agnitu des) and J (Vega -based magnitudes) are from the quasar discovery papers 
Fan et al.l (|200ol ). iFan et al.l (|200ll ). and lFan et all (|2004l ). 



-29- 



Table 2. Redshifts and Continuum Slopes 



Quasar (SDSS) Redshift Slope uv (a I/ ) a Slope pt(a^) b 



J0836+0054 


5.810 


± 


0.003 c 


-0 62+ - 06 
u - DZ -0.06 


-0.52 


J1030+0524 


6.309 


± 


0.009 


46 +a18 


-0.31 


J1044-0125 


5.778 


± 


0.005 d 


-n 27+ ' 09 

u - z ' -0.10 


-0.30 


J1306+0356 


6.016 


± 


0.005 


50+ ' 12 
u.ou_ 14 


-0.12 


J1411+1217 


5.927 


± 


0.004 


-0 21 +a16 

u - zi -0.19 


-0.51 


J1623+3112 


6.247 


± 


0.005 




-0.32 



a UV continuum slopes measured from the Gemini spectra. 
b Optical continuum slopes taken from I Jiang et al.l (2006). 



c Taken from lKurk et al.1 ([20071). 

d Measured from the Cuil A1909 emission line. 



Table 3. Emission Line Properties 



J0836+0054 J1030+0524 J1044-0125 J1306+0356 J1411+1217 J1623+3112 



Nv 


EW 


15.9±0.7 




FWHM 


31.5±1.0 


Sin 


EW 


3.6±0.3 




FWHM 


15.0±0.6 


Oi 


EW 


7.4±0.1 




FWHM 


24.6±0.5 


Cn 


EW 


5.8±0.5 




FWHM 


41.9±3.1 


Si iv 


EW 


9.3±4.2 




FWHM 


25.3±5.9 


C, TV 

\J 1 V 


EW 

J_J vv 


i rc n_i_i f. 

_L O .\J— 1 — L .U 




FWHM 


37.1±3.0 


Hen 


EW 






FWHM 




Om] 


EW 






FWHM 




Aim 


EW 


1.8±0.9 




FWHM 


16.3±6.8 


Cm] 


EW 


18.9±2.3 




FWHM 


53.2±4.9 


Mgn 


EW 






FWHM 





4.9±0.1 3.9±0.2 

8.7±0.2 10.4±0.3 

3.0±0.2 2.3±0.1 

15.8±0.7 12.7±0.3 

6.7±0.3 2.6±0.2 

18.8±0.7 24.3±1.3 
4.0±0.4 
17.0±1.5 

8.5±1.7 8.9±0.7 

18.4±2.9 35.7±1.8 

40.0±3.1 24.7±2.1 

27.3±2.9 42.2±3.1 

3.4±2.0 

21.7±8.8 

2.3±1.3 

22.9±11.0 

5.0±1.6 
31.0±6.7 

19.3±2.8 24.0±1.8 

31.6±3.3 50.3±3.2 

46.1±3.7 

33.8±2.2 



9.3±0.3 28.9±0.4 
15.0±0.2 22.3±0.3 
1.7±0.1 0.9±0.1 
12.8±0.4 3.7±0.3 
6.7±0.2 
8.1±0.2 



9.5±0.8 8.1±6.4 

29.0±2.0 21.5±12.2 

29.6±2.4 76.2±5.8 

24.5±2.3 17.5±1.6 

2.6±1.8 5.1±2.5 

15.5±7.5 39.0±26.8 

2.7±2.0 2.8±2.5 

26.9±19.4 42.2±34.0 

7.0±3.1 5.1±2.3 

34.5±14.7 25.6±9.2 

19.1±3.9 26.1±3.2 

37.6±7.4 31.3±2.6 

24.9±2.1 37.5±5.5 27.4±3.9 

33.2±2.2 23.7±2.7 34.4±2.4 



Note. - 



Rest-frame FWHM and EW are in units of A. 



Table 4. Flux Ratios 



J0836+0054 J1030+0524 J1044-0125 J1306+0356 J1411+1217 J1623+3112 



Nv/Civ 

Sin/Civ 

Oi/Civ 

Cii/Civ 

Siiv/Civ 

Hcii/Civ 

O m] /Civ 

Aim/Civ 

Cm] /Civ 

Feii/Mgn 



1.44±0.17 
0.32±0.04 
0.62±0.07 
0.47±0.06 
0.72±0.33 



0.10±0.05 
0.95±0.15 



0.21±0.02 
0.12±0.02 
0.26±0.02 
0.14±0.02 
0.27±0.06 
0.07±0.04 
0.05±0.03 

0.29±0.05 
5.45±0.58 



0.23±0.03 
0.13±0.02 
0.14±0.02 

0.43±0.05 



0.15±0.05 
0.68±0.08 



0.55±0.05 
0.10±0.01 



0.41±0.05 
0.08±0.05 
0.08±0.06 
0.15±0.07 
0.38±0.08 
5.77±0.72 



0.56±0.04 
0.02±0.01 
0.12±0.02 

0.13±0.10 
0.06±0.05 
0.03±0.02 
0.05±0.02 
0.24±0.03 

5.27±0.72 3.21±0.62 
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Table 5. Central BH Masses (10 9 M ) 



Quasar (SDSS) 


LBol d 


LBol/LEdd b 


Mbh(Civ) 


M BH (Mgn) 


M BH (Mgn) c 


J0836+0054 


47.72 


0.44 


9.3±1.6 






J1030+0524 


47.37 


0.50 


3.6±0.9 


1.0±0.2 


2.1±0.4 


J1044-0125 


47.63 


0.31 


10.5±1.6 






J1306+0356 


47.40 


0.61 


3.2±0.6 


l.lrbO.l 


2.2±0.3 


J1411+1217 


47.20 


0.94 


1.3±0.3 


0.6±0.1 


0.9±0.2 


J1623+3112 


47.33 


1.11 




1.5±0.3 





a Bolometric luminosity in log[erg s x ] from I Jiang et al.1 (2006). 

b LEdd is derived from Mbh(Civ) except for SDSS J1623+3112, whose L^dd is 
derived from MBH(Mgn). 

c Mg H is estimated from the new relation by Vestergaard et al. (in preparation). 
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Table 6. Absorption Lines 



Absorption system 


Kbs (A) 


ID 




W (A) 


JUooO-|-UU04 


loZDl.ol 


Mb- tt 


(2796) 


q 7/1 OR 


2. Oil 




lozyt>.4o 


Mgll 


(2803) 


6. 14:24: 


9 nc 
z.Uo 


SDSS J1306+0356 (a) 


8952.03 


Mgn 


(2796) 


2.2013 


1.86 




8977.30 


Mgn 


(2803) 


2.2022 


1.99 




9137.13 


Mgi 


(2852) 


2.2028 


1.29 


SDSS J1306+0356 (b) 


9877.64 


Mgn 


(2796) 


2.5323 


3.24 




9903.27 


Mgn 


(2803) 


2.5324 


3.15 




9185.23 


Fen 


(2600) 


2.5325 


1.94 


SDSS J1306+0356 (c) 


16405.55 


Mgll 


(2796) 


4.8668 


2.92 




16447.67 


Mgn 


(2803) 








9795.84 


Alii 


(1670) 


4.8632 


1.11 




15252.21 


Fen 


(2600) 


4.8658 


1.49 


SDSS J1306+0356 (d) 


16448.85 


Mgn 


(2796) 








16491.08 


Mgll 


(2803) 


4.8823 


2.11 




9822.93 


Alii 


(1670) 


4.8795 


0.76 




15293.94 


Fen 


(2600) 


4.8819 


1.97 



